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Abstract 

We consider the semiclassical quantization of the Sinai billiard for disk 
radii R small compared to the wave length 2-Tr/k. Via the application of 
the periodic orbit theory of diffraction we derive the semiclassical spectral 
determinant. The limitations of the derived determinant are studied by 
comparing it to the exact KKR determinant, which we generalize here for the 
Ai subspace. With the help of the Ewald resummation method developed for 
the full KKR determinant we transfer the complex diffractive determinant to 
a real form. The real zeros of the determinant are the quantum eigenvalues in 
semiclassical approximation. The essential parameter is the strength of the 
scatterer c = Jo(kR)/Yo(kR). Surprisingly, this can take any value between 
— oo and +00 within the range of validity of the diffractive approximation 
kR <C 4, causing strong perturbation in spite of the smallness of the disk. We 
study the statistics exhibited by spectra for fixed values of c. It is Poissonian 
for c = ±00, provided the disk is placed inside a rectangle whose sides obeys 
some constraints. For c = we find a good agreement of the level spacing 
distribution with GOE, whereas the form factor and two-point correlation 
function are similar but exhibit larger deviations. By varying the parameter 
c from to ±00 the level statistics interpolates smoothly between these 
limiting cases. Any transitional level statistics can thus be found in the 
spectrum of a Sinai billiard with sufficiently small R while we go from the 
quantum to the semiclassical limit k — » 00. 
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1 Introduction 



One of the fundamental questions of quantum chaos and mesoscopic physics is 
the relation between the classical motion and the energy levels of the quantum 
counterpart of the system |1|, 0, |3|. Correlation of energy levels in chaotic systems 
in the semi classical limit (% — > 0) are expected to coincide with those in Gaussian 
Random Matrix Ensembles (RME) with different symmetries (GOE, GUE etc.), 
while generic integrable systems are characterized by Poissonian level statistics ||12||. 
Recently significant progress has been made concerning the applicability of RME for 
disordered and chaotic systems^. In disordered systems the two point correlation 
function has been expressed in terms of the spectral determinant of the diffusion 
equation. These results have been derived with supersymmetric techniques. It has 
been suggested that this result can be generalized to chaotic systems ||, [|, where 
the spectral determinant of the Perron- Frobenius operator plays a similar role. 

Similar results have been derived with the help of semiclassical Gutzwiller type 
periodic orbit expressions^. This approach relies heavily on the so-called diagonal 
approximation for t <t*, where t* is of the order of the Heisenberg time (see below), 
in order to, in an ingenious way, make statements about t > t*. The suggested 
results of 0, |5[] then follows after some further approximations. 

One of the conclusions of these studies is that deviations from RME behavior is 
expected only when the chaotic mixing is slow and the characteristic time of mixing 
T m is comparable with the Heisenberg time. In a usual system the mixing time T m 
is a classical quantity characterizing the decay of classical correlation functions and 
as such it is of order h°. The Heisenberg time is the uncertainty time Th ~ h/A 
corresponding to the mean level spacing (A) . In the semiclassical limit h — > the 
mean level spacing behaves like A ~ h d } where d is the dimension of the system 
and the Heisenberg time gets much larger than the mixing time. 

Slow mixing can be expected if we perturb an integrable system with a small 
point like object, whose classical size can be neglected. In this case the only source 
of mixing is the quantum mechanical diffraction of waves on the small object. The 
diffraction constant governing the speed of the mixing is proportional with h 2 and 
the resulting mixing time is proportional with 1/h. In a two dimensional system 
the level spacing is also proportional with h 2 and the ratio of the mixing time and 
the Heisenberg time remains fixed in the semiclassical limit 7i — > 0. This gives 
the possibility to study systems, where corrections to RME remain finite. In this 
paper we study the Sinai billiard with a small scatterer, which, as we show here, 
is a strongly diffractive system and its level statistics deviates strongly from the 
RME result. 

We are considering the standard Sinai billiard, a unit square with a reflective cir- 
cular disk in the middle. We then develop the periodic orbit description of the Sinai 
billiard with small scatterer based on the diffractive extension |7]] of the Gutzwiller 
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trace formula. This is done in sections 2 and 3.1. The diffractive determinant we 
derive from the trace formula is directly related to the KKR determinant of Berry, 
which is discussed in section 3.2. We will study a straightforward generalization of 
Berry's method to the Ai subspace, which exhibit some very interesting features, 
even in the diffractive or small radius limit. In section 4 we study the level statis- 
tics exhibited by the zeroes of the diffractive determinant for various values of the 
parameter kR. We then change from a square to a suitably chosen rectangle shape, 
in order to avoid the non generic degeneracies of the square billiard. 

2 Geometric theory of diffraction 

In this section we review the main steps of the derivation of the diffractive trace 
formula and spectral determinant introduced in Ref . [|7| . The emphasis is on point- 
like diffractive sources with no angular dependence in their associated diffraction 
constants. These point sources are placed inside a bounded and otherwise non 
diffractive system. 

2.1 Green functions 

For the description of diffraction effects in the semiclassical limit (h — > 0) Keller § 
introduced the Geometric Theory of Diffraction (GTD). According to this theory, 
the energy domain Green function can be divided into a geometric (G) and a 
diffractive (D) part: 

G(q,q',E) = G G (q,q',E)+G D (q,q',E) , (1) 

where the geometric part is the semiclassical Green function. The semiclassical 
Green function is a sum for all classical trajectories (j) with energy E connecting 
the starting point q and the end point q' 

G G (q,q',E)= J2 G { j\q,q',E) , (2) 

jeq-^q' 

where Go^(g, q', E) is the Van Vleck-Gutzwiller Green function. It can be calculated 
from the action Sj(q, q', E) = / 9 pdq calculated along trajectory j 

GH\q,qi,E)=A 3 (q,q',E)eT^'> E ) , (3) 

where the amplitude Aj(q,q' , E) is the Van Vleck-Gutzwiller determinant f| mul- 
tiplied with the Maslov phase. 

The diffractive part of the semiclassical Green function Go(q,q',E) describes 
the effect of wave scattering on points qf. being singular points of the classical 
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dynamics^, PH . It is a sum over all possible trajectories which start from q, go 
to one of the singular points q s k1 then leave the singular point and reach the end 
point q'; or go and visit any combination of singular points before reaching the end 
point. A trajectory going trough one singularity contributes to the Green function 
withg 

Gi 3 \q,q!,E)d k (p,p>,E)G ( /\q!,q',E) . (4) 

The generalization to multi scattered trajectories is straightforward, and will be 
developed in the next section. In eq (|j), j and j' are indices of the incoming and 
the outgoing classical trajectory, d k (p,p',E) is the diffraction constant associated 
with the singularity located at qf.. It depends on the direction of incoming and 
outgoing momenta p and p' and energy E. The diffraction constant is the leading 
part of the exact quantum mechanical scattering amplitude of the singularity. 



2.2 Traces and periodic orbits 

Based on the expression for the Green functions in the previous subsection we can 
derive the diffractive version of the Gutzwiller trace formula]^]]. 

The trace of the total Green function is a sum of its geometric and diffractive 
parts: 

TrG(E) = TrG G {E) +TrG D {E) + g (E) , (5) 

where the geometric part TrGo(E) is just the usual Gutzwiller trace formula, a 
sum over unstable periodic orbits possibly amended by a sum over neutral orbits 
and the contribution from zero length orbits is in go(E). Gg{0, 0, E) is essentially 
the Green function Gg{0, 0, E) without the singular contribution coming from the 
zero length orbit. We concentrate here on the new term 

TrG D (E) = jG D (q,q,E)dq . (6) 

To be more specific, we assume that the diffraction is caused by the presence 
of N small disks located at g^, where 1 < k < N, whose diffraction constants 
dk(E) do not depend on the scattering angle [11]. This will simplify the formulation 



considerably. We can easily see that the contribution to the Green function from 
all orbits via the disks k\ and 

E Gf\q,q s kx )d kx J2 G^\q s kv qt 2 )d k2 ... £ G^\q^q') (7) 

= E G^\q,q s ki )d kl G G (q s ki ,q s k2 )d k2 ... £ G^\q^q') . 
The generalization to multi scattered trajectories is now obvious. 
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We can introduce a symbolic dynamics by simply labeling the disks by the cor- 



responding indices k. In an open N-disk system a periodic disk sequence k\k 2 ■ ■ ■ k n 
codes the single periodic orbit visiting the disks in the prescribed order. In a 
bounded system it denotes all such periodic orbits e.g. the transition . . . kik i+ i . . . 
denotes all (topologically distinct) ways to go from disk ki to disk k i+ i. In the fol- 
lowing we will refer to periodic symbol codes and symbolic periodic orbits to avoid 
confusion with periodic orbits. 

The calculation of the trace of the Green function Gjj{q, q', E) amounts to add 
up contributions from all distinct paths from q to q' via n > 1 disk scatterings 
and then taking the trace. The basic building block Gc{ki, fci+i, E) in semiclassical 
approximation is a sum over all paths from (see (0)). 

Consider now a situation where the starting and ending points coincide q = q'. 
The contribution from the symbolic path k\ , k 2 . . . k n to the trace is 

fdq £ G^\q,qtJd kl G G (ql i ,q^d k2 ---G G (ql^ql)d kn £ G^\q s kn ,q) 

(8) 

The integral over q can be performed in the H — > limit by using the saddle point 
method. The result is 



•^G(g^,^)4iG G (^,^)d^---G G (g^_ 1> gJ n ) d *» > ( 9 ) 

The stationary phase conditions select out all points q lying on a trajectory between 
disk k n and k±, that is, it selects periodic orbits. 

To get the full contribution from the symbolic periodic orbit k\k 2 ■ ■ -k n one must 
add the contribution from all nonidentical cyclic permutations of the sequence 
k\k<i . . . k n in eqs. (|8|) and (|Sp. The symbolic periodic orbit k\ki ■ ■ ■ k n can be a 
repetition of a shorter primitive (symbolic) periodic orbit k\k 2 . . . k n = jf. The 
length of the code p is denoted n p = n/r. Instead of restricting ourself to all non 
identical permutations we can take all cyclic permutations and then divide by r. 

The result of this somewhat laborious exercise in combinatorics can be conve- 
niently expressed in terms of the derivative of the quantity 



t p = II dkpciqk^qk^) (10) 

3=1 



as 

00 1 d 



P £fL D r=l ' arj 
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where Qn is the set of all primitive periodic symbol sequences taken from the 
alphabet {k; 1 < k < N} where N is the number of singular points in the system. 

We now introduce the spectral determinant A(E), defined through the loga- 
rithmic derivative 

TrG = -^\ogA(E) . (12) 

Since the Green function is a sum (see Eq.|5|), the spectral determinant naturally 
factorizes product 

A(E) = A (E) ■ A G (E) ■ A D (E) , (13) 

where A (E) is associated with zero length orbits and corresponds to the mean 
level density, Ag{E) is the Gutzwiller-Voros zeta function associated with classical 
periodic orbits, possibly amended by neutral orbits. The diffractive determinant 
(or zeta function) now reads 

A D (E)= J] (1-fp) • (14) 

3 Application to the Sinai billiard 

In this section we work out the diffractive spectral determinant of the Sinai billiard 
in the small disk limit. Then we show that the diffractive determinant is essentially 
the leading part of the Korringa-Kohn-Rostoker (KKR) determinant in the small 
disk limit. 



3.1 The diffractive spectral determinant 



d(E) ~ -4*7^7 • (15) 



The application for the Sinai billiard is now very simple, provided kR is sufficiently 
small. The diffraction constant has been derived in Ref. ||11|| , in the limit kR — > 
it was found to be 

Jo(kR) 
Hg-\kR) 

It does not depend on the opening angle and the considerations in the previous 
section can be applied. 

The billiard contains only a single disk and the symbolic dynamics introduced 
in the previous section has only a single letter implying the existence of one single 
primitive periodic code, and the resulting determinant is trivial 

A D (E) = l-d(E)G G (0,0,E) . (16) 

The geometry of the Sinai billiard is easy to visualize if one unfolds it into a 
quadratic lattice of disks. A trajectory from the diffractive object (i.e. the small 
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disk) and back to itself is characterized by a lattice vector p = (p x , p y ) where p x 
and p y can take any integer value. The geometric Green function Gq(0, 0, E) may 
be expressed as a sum over lattice vectors 

G G (0,0,E) = ]T G (0,p,E) , (17) 

p^(0,0) 

with the zero length vector p = (0, 0) excluded. The Green function Go may be 
expressed in terms of (the semiclassical limit of) a free Green function 

G G (0,0,£) = £ G free (0,P,E)(-l) p * +p « , (18) 

p^(0,0) 

where (— iy x+ Pv is a phase index associated with the reflection on the square wall. 
But in order to improve the approximation we use the exact free Green function 
in two dimensions 

G free (x,x') = - l -H^\k\x-x'\) . (19) 
The resulting expression is 

A D (E) = l-d(E) £ {- l -)H$\k\p\){-iy^ , (20) 

^4(0,0) 

where k = \JlE is the momentum in dimensionless units h = m = 1. The sum for 
the lattice vectors is the exact Green function of the empty square with the zero 
length term removed. 

The Green function can also be expressed via the eigenfunctions of the empty 
square. The wave functions can be categorized according to the different irreducible 
representations of C± v . There are four one dimensional (Ai, A2, B\ and B2) repre- 
sentations and one two dimensional (E). The disk is centered at the point where 
all the symmetry lines meet. Wave functions of the empty square billiard vanish 
at this point except for the A\ subspace. Accordingly, only the energies belonging 
to this subspace are affected by diffraction in leading semiclassical order as will be 
demonstrated later. 



3.2 The KKR determinant 

The full quantum Sinai billiard problem can be solved with the KKR method 
introduced in this context by Berry in Ref. [|12[|, where the A 2 subspace has been 
worked out in detail. The method leads to a determinant D^kr whose zeroes 
yield the eigenvalues of the problem. The generalization to other subspaces are 
fairly straightforward. Here, without repeating the derivation of Ref.[12|], we quote 
the relevant results for the A\ and A2 subspaces and show that the periodic orbit 
expression (BO) is essentially the leading part of the KKR determinant. 



7 



The effective dimension of the KKR matrix is ~ (kR)/4. For kR < 4, which 
we call the diffractive zone, it is essentially one dimensional. The determinant is 
given by the leading element of the matrix 

a fu\ - J 1 + 48§S r (e) in the A 1 subspace 

KKR{ ) " \ 1 + gjgf(S5(e) - SJ(e)) in the A 2 subspace ' {21) 

where the real functions Sj r ^(e) are related to the structure constants 



via 



S , ,(e) = 5r(e)+i<Jo,i • (23) 

e = (/c/27r) 2 is a rescaled energy (denoted by E in fl2fl). 

The structure constants, as they stand, do not converge. They can be evaluated 
by transforming them to sums over the dual lattice using the Ewald resummation 
method [P^j. The resulting expressions are 

s j %)K ± u ^ eMm-f/e)) ( 
r7 e e — v z 

n z ^7 e — v z 
We have introduced the vector v in the dual lattice 



v = (v x + |, Uj, + |) in the Ai subspace 
v = (i> z , t> y ) in the A 2 subspace 



(26) 



where v x and f y are integers. The poles of the expressions located in v 2 are the 
eigenvalues of the empty square billiard. The smaller the number Q the more 
accurate the method works. 



Combining Eqs. (|15j ), (|20f) , (f21|) and (|22|) we find the desired relation between 
the diffractive limit of the KKR determinant and the diffractive determinant 



^'I'-Sj^ 1 ' (27) 

This relation is one of the main results of the paper. 

The KKR determinant, and thus the diffractive spectral determinant, exhibit 
poles at the unperturbed eigenvalues e^ -* = v 2 . They cancel the zeroes of the 
geometric determinant Ac(fe) according to the formula 

A(k) ~ Aa(k)A D (k) , (28) 
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and give rise to corrected eigenvalues given by the zeroes of Ai>(k). This is an 
important property of diffractive zeta functions: they cancel the geometric zeros 
and create new ones. 

Note that there is no contribution coming from isolated unstable orbits in our 



system, log Ac(k) is just the Berry- Tabor sum for the square billiard [13]. 



4 Investigation of the energy levels 

The investigation of the KKR determinant for all symmetry subspaces is not a 
superfluous exercise. Next we show that the Ai subspace will exhibit some very 
interesting features, to which we devote the rest of the paper. 

The quantum eigenvalues in the diffractive limit can be determined from the 
equation 

A D (fc) = 0. (29) 
or euivalently from the real equation 

S (r) {e) + c{kR) =0 , (30) 

where 

c(kR) = ^§ • (31, 
J (kR) 

4.1 Perturbation theory 

When c(kR) is large one can use perturbation theory to compute the roots of Eq. 
(|30| ) (see Ref . [|T2f ) . One thus assumes that the difference between a perturbed 
level e and an unperturbed level is small. Sq (e) close to may then be 
approximated by 



7r^ e 



yielding the perturbed eigenvalue 



eKie(0) _l^m (33) 

tH Y (kR) 

Fig. 1 illustrates how well this expression works for R = 0.025, where we have 
plotted the distance from an eigenvalue of the Sinai billiard (R = 0.025) and the 
nearest eigenvalue of the square billiard. The expression above applies only to 
non-degenerate states. For a n-degenerate unperturbed energy, n — 1 levels remain 
almost unaffected, which is clearly seen on Fig. 1. 

There are two regions in the range < kR < 4 where perturbation theory 
should work in principle. 
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The first is the limit kR — > 0. Taking well known asymptotic expressions [H 
for the Bessel functions one gets 

e = e (0) - - - (34) 

2 7 rlog(A; J R/2)+7 ' 1 ' 

where 7 = 0.57721 ... is Euler's constant. However, this expression goes very slowly 
to zero as kR — > and unless the radius is very small perturbation theory does not 
even apply to the ground state = 0.5. 

The second is where c(kR) has a pole within the diffractive zone kR < 4, namely 
at kR = jo,i = 2.40482 . . .. Close to this value, the perturbation expression works 
very well, as it can be seen on the scatter plot Fig. 1. Note, that the index of 
the levels is now shifted one step up relative to the unperturbed ones, this is why 
we chose to plot the nearest distance to an unperturbed eigenvalue. We note that 
c{kR) exhibit a zero at kR = y ,i ~ 0.899 so the perturbation theory also breaks 
down deep inside the diffractive zone. 



As a comparison, in the A 2 subspace the perturbation is on average |fL2 




8 J t (kR) , ... 
<e - e ° > =-^YjkRj~-^{-) ■ <35) 

Perturbation theory works well almost throughout the whole diffractive region. 
So the spectrum in the A2 representation, in the diffractive zone, can hardly be 
expected to be particularly interesting. 

For A\, on the other hand, \c(kR)\ will range from to 00 in the interval 
< kR < ?/o,i ~ 0.899 which is well inside the diffractive zone. 



4.2 Level statistics in the diffractive limit 

In this section we solve Eq. ( |3"0"| ) numerically and study the statistics of levels. 
The essential parameter of the problem is kR since the mean level spacing A e does 
not depend on energy e. Here we will study spectra for fixed values of c(kR). 
The idea of keeping c and thus kR fixed needs some explanation. If we decrease 
R and keep kR fixed by increasing k, more and more levels are squeezed in a 
neighborhood of some k = k^. We then can safely replace kR by koR and it is 
possible to obtain a sufficiently large number of levels around ko for which the 
value of c is approximately constant and given by c(koR). For obvious reasons we 
calculate the bottom part of the spectrum. However, this sample should exhibit the 
same statistics as any comparable sequence of levels, e.g. around k . The reason 
is that the location of a level depends on the parameter kR and the location of 
the unperturbed levels in the environment and the statistics of these unperturbed 
levels is well defined and given by a Poissonian process. 
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In the actual calculation we made a minor modification of the conventional Sinai 
billiard, since the empty square billiard has degeneracies in its spectrum. Their 
origin is number theoretical and is connected with the high degree of symmetry in 
the system. If a disk is present, these degeneracies are split by higher order terms 
in the KKR determinant. Within the scope of the present diffractive determinant 
they remain there. In the following we want to avoid this complication and replace 
the square by a rectangle with sides a = 1 and b = 2 1 / 4 . The function Sq (e) is 
modified accordingly 

x 1 1 ^ exp(g[l - v 2 /e]) 

S o(e)^——}^ Ei(Q)ir, 36 

ao ir z e — v z 

where v = ((v x + 1/2) /a, (v y + l/2)/b) is now the dual lattice of the rectangle. 

We investigate three statistical measures on the spectra. First we will investi- 
gate the integrated probability distribution of neighboring levels 

P( s ) = J26(s-(e n+1 -e n )/A e ) , (37) 

n 

where A e is the mean level spacing around the energy e and 9(x) is the step 
function. Then we investigate the two point correlation function of levels 

R(e) = ( £ <K( e - e t )/A e + e/2) ■ 5((e - e,)/A e - e/2) ) e (38) 

ij 

= 5(e) + ( £ 8({e - e,)/A e + e/2) • 8{(e - e 3 )/A e - e/2) ) e = 5(e) + R(e) 

where the average is taken for a large number of energies. We also compute its 
Fourier transform, the form factor, 

/+oo 
e 2mrt (R(e) - l)de (39) 

The correlation function and the form factor are computed over a gaussian window 
centered at the middle of the sample spectrum, its width is about one sixth of the 
sample size. The results are then smeared with another gaussian. We choose the 
width of this gaussian rather small so some (non significant) oscillations due to 
finite statistics will remain in the plots. 

In Fig. 2 we show results for c = 0, with 5000 levels computed. The level spacing 
curve (Fig 2a) is close to the prediction of the random matrix theory (GOE), with 
a small but significant deviation. An exact agreement is hardly possible because 



the levels are locked between the unperturbed levels|12|. The close agreement with 
GOE is in this perspective rather surprising. 

The formfactor K(t) is plotted in Fig 2b. There is again good agreement with 
GOE except for small r. The bump for small values of r indeed persist rather 
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high up in the spectrum (of the full KKR determinant) and is quite similar to the 
behavior predicted in [p~5j| . 

The correlation function R(e) is plotted in Fig. 2c. The deviation from GOE 
is now most pronounced around e ~ 1. 

In Fig. 3 we show the correlation function -R(e) for a geometric sequence of c's, 
each spectrum contains 1500 levels. We note how it, in a uniform way, interpolates 
between a Poissonian result and the GOE-like. 



5 Concluding Remarks 

The KKR determinant may be recast into a periodic orbit expressions in the limit 
kR — > oo |I2|| . Periodic orbit calculations are very laborious in this limit, it is 
questionable whether it provides a useful method of extracting (many) eigenvalues 
for bound chaotic systems |T6[ . 

What motivated us in this paper, is that the KKR determinant yields (diffrac- 
tive) periodic orbit expressions also in the limit kR — > 0. This diffractive limit is 
very feasible form a periodic orbit point of view. 

A traditional symbolic dynamics for the Sinai billiard uses the lattice vectors p 
as symbols [17]]. In the diffractive limit, this symbolic dynamics is complete. The 



reason for this is the absence of shadow behind disks in the kR — > limit. Moreover, 
the diffraction constant is independent of the scattering angle. In the language of 
cycle expansions ||18||, these circumstances imply that the cycle expansion (|20"D is 
entirely given by the (well defined) fundamental part, the curvature correction 
vanish identically. 

So, what we have presented here is a very straight forward application of the 
cycle expansion technique as well as of the geometric theory of diffraction. Yet, 
the results are non-trivial. Even if we have not used the diffractive determinant 
for computations as it stands ( p0[) we have demonstrated that study of spectral 
statistics is in principle approachable for periodic orbit theories in this case. 

A natural extension of this work is to increase the number of small scatterers 
and study whether GOE emerges as the number of scatterers tends to infinity and 
are distributed randomly, preliminary results indicate that this is indeed the result 
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Figure 1: Distance between eigenvalues of the Sinai billiard and the nearest eign- 
value of the empty square compared with perturbation theory. The eigenvalues are 
from the A 1 subspace for radius R = 0.025. 
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Figure 2: a) Integrated level spacing distribution, b) Form factor, c) Two point 
correlation function. Eigenvalues computed from eq ( p0|) with c = for a rectan- 
gular Sinai billiard. 
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